Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Role of Equatorial Anomaly in Earthquake time precursive features: A few strong events over West Pacific zone

The earthquake (EQ) time coupling processes between equator-low-mid latitude ionosphere are complex due to inherent dynamical status of each latitudinal zone and qualified geomagnetic roles working in the system. In an attempt to identify such process, the paper presents temporal and latitudinal variations of ionization density (foF₂) covering 45°N to 35°S, during a number of earthquake events (M?>?5.5). The approaches adopted for extraction of features by the earthquake induced preparatory processes are discussed in the paper through identification of parameters like the ‘EQ time modification in density gradient’ defined by δ?=?(foF_{2 max}???foF_{2 min})∕τ_mm, where τ_mm – time span (in days) between EQ modified density maximum and minimum, and the Earthquake time Equatorial Anomaly, i.e. EEA, one of the most significant phenomenon which develops even during night time irrespective of epicenter position. Based on the observations, the paper presents the seismic time coupling dynamics through anomaly like manifestations between equator, low and mid latitude ionosphere bringing in the global Total Electron Content (TEC) features as supporting indices.     

Efficiency Estimation of Electrothermal Thrusters Use in the Control System of the Technological Spacecraft Motion

The influence of various control systems of the orbital motion of a technological spacecraft on the level of microacceleration of its internal environment is simulated. Conclusions are drawn about the effectiveness of control systems with different actuators for realization of certain gravitationally sensitive processes onboard a spacecraft.     

Diurnal and seasonal variations in long-term changes in the E-layer critical frequency

Long-term changes in the E-layer critical frequency, foE, at three stations of the European region (Juliusruh, Slough and Rome) and also at Moscow and Wakkanai stations are analyzed by the method developed by the authors and described in detail in the previous papers. It is found that Juliusruh and Slough stations demonstrate a well-pronounced change in foE (a trend) during two previous decades. At the same time, the same features of the behavior of the aforementioned trend k(foE) are obtained. The trend is positive and negative in the morning and evening hours, respectively. Similar diurnal behavior of k(foE) is found also for Moscow station but with lower absolute values of the trends. A well-pronounced seasonal behavior of k(foE) is detected at Juliusruh and Slough: the trend is minimal and maximal in the summer period and at the end of fall—beginning of winter, respectively. The maximal amplitude in the morning hours reaches +0.04?MHz per year, whereas the minimal amplitude in evening hours is ?0.06?MHz per year. No systematic changes exceeding by the magnitude 0.01?MHz per year are found for Rome and Wakkanai stations. It is assumed that the observed trends are related to changes (trends) in the meridional wind bringing NO molecules from the auroral oval to lower latitudes.     

Assessment study of ionosphere correction model using single- and multi-shell algorithms approach over sub-Saharan African region

It is a known fact that ionosphere is the largest and the least predictable among the sources of error limiting the reliability and accuracy of Global Navigation Satellite Systems (GNSS) and its regional augmentation systems like Satellite Based Augmentation System (SBAS) in a safety-of-life application. The situation becomes worse in the Equatorial Ionization Anomaly (EIA) region, where the daytime ionization distribution is modified by the fountain effect that develops a crest of electron density at around ±15° to ±20° of the magnetic equator and a trough at the magnetic equator during the local noon hours. Related to this phenomenon is the appearance of ionosphere irregularities and plasma bubbles after local sunset. These may degrade further the quality of service obtained from the GNSS/SBAS system of the said periods. Considering the present operational augmentation systems, the accuracy and integrity of the ionosphere corrections estimate decreases as the level of disturbances increases. In order to provide a correct ionosphere correction to the user of GNSS operating in African EIA region and meet the integrity requirements, a certified ionosphere correction model that accurately characterizes EIA gradient with the full capacity to over-bound the residual error will be needed. An irregularities detector and a decorrelation adaptor are essential in an algorithm usable for African sub-Saharan SBAS operation. The algorithm should be able to cater to the equatorial plasma vertical drifts, diurnal and seasonal variability of the ionosphere electron density and also should take into account the large spatial and temporal gradients in the region. This study presents the assessment of the ionosphere threat model with single and multi-layer algorithm, using modified planar fit and Kriging approaches.     

Precise line-of-sight modelling for angles-only relative navigation

This work presents a precise analytical model to reconstruct the line-of-sight vector to a target satellite over time, as required by angles-only relative navigation systems for application to rendezvous missions. The model includes the effects of the geopotential, featuring: the analytical propagation in the mean relative orbital elements (up to second-order expansion), the analytical two-way osculating/mean orbital elements’ conversion (second-order in $J_2$ and up to a given degree and order of the geopotential), and a second-order mapping from the perturbed osculating elements’ set to the local orbital frame. Performances are assessed against the line-of-sight reconstructed out of the precise GPS-based positioning products of the PRISMA mission. The line-of-sight modelled over a far-range one day long scenario can be fitted against the true one presenting residuals of the order of ten arc-seconds, which is below the typical sensor noise at far-range.